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Abstract 

A game theoretic approach is used to derive the optimal decentralized power allocation (PA) in fast fading 
multiple access channels where the transmitters and receiver are equipped with multiple antennas. The players 
' (the mobile terminals) are free to choose their PA in order to maximize their individual transmission rates (in 

l/^ , particular they can ignore some specified centralized policies). A simple coordination mechanism between users 

in 

i is introduced. The nature and influence of this mechanism is studied in detail. The coordmation signal indicates 

Csl ' to the users the order in which the receiver applies successive interference cancellation and the frequency at 

O ' 

I which this order is used. Two different games are investigated: the users can either adapt their temporal PA to 



o 

> 



their decoding rank at the receiver or optimize their spatial PA between their transmit antennas. For both games 
a thorough analysis of the existence, uniqueness and sum-rate efficiency of the network Nash equilibrium is 



X 

^ • conducted. Analytical and simulation results are provided to assess the gap between the decentralized network 



performance and its equivalent virtual multiple input multiple output system, which is shown to be zero in some 
cases and relatively small in general. 
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I. Introduction 

We consider a special case of decentralized or distributed wireless networks, the decentralized multiple access 
channel (MAC). In this context, the MAC consists of a network of several mobile stations (MS) and one base 
station (BS). In the present work, the network is said to be decentralized in the sense that each user can choose 
freely his power allocation (PA) policy in order to selfishly maximize a certain individual performance criterion. 
This means that, even if the the BS broadcasts some specified policies, every (possibly cognitive) user is free 
to ignore the policy intended for him if the latter does not maximize his performance criterion. 

The problem of decentralized PA in wireless networks is not new and has been properly formalized for the 
first time in ||2l, ||3l. Interestingly, this problem can be formulated quite naturally as a non-cooperative game 
with different performance criteria (utilities) such as the carrier-to-interference ratio ||4l, aggregate throughput 
|[5l or energy efficiency 161, n. In this paper, we assume that the users want to maximize information-theoretic 
utilities and more precisely their Shannon transmission rates. Many reasons why this kind of utilities is often 
considered are provided in the literature related to the problem under investigation (some references are provided 
further). Here we will just mention three of them. First, Shannon transmission rates allow one to characterize 
the performance limits of a communication system and study the behavior of (selfish) users in a network where 
good coding schemes are implemented. As there is a direct relationship between the achievable transmission 
rate of a user and his signal-to-interference plus noise ratio (SINR), they also allow one to optimize performance 
metrics like the SINR or related quantities of the same type (e.g., the carrier-to-interference ratio). From the 
mathematical point of view. Shannon rates have many desirable properties (e.g., concavity properties), which 
allows one to conduct deep performance analyses. Therefore they provide useful insights and concepts that are 
exploitable for a practical design of decentralized networks. Indeed, the point of view adopted here is close to the 
one proposed by the authors of Il8l for DSL (digital subscriber lines) systems, which are modeled as a parallel 
interference channel; lH for the single input single output (SISO) and single input multiple output (SIMO) 
fast fading MACs with global CSIR and global CSIT (Channel State Information at the Receiver/Transmitters); 
IfTOl for MIMO (Multiple Input Multiple Output) MACs with global CSIR, channel distribution information at 
the transmitters (global CDIT) and single-user decoding (SUD) at the receivers; ifTTI . |[T2l for Gaussian MIMO 
interference channels with global CSIR and local CSIT and, by definition of the conventional interference 
channel |[T3l . SUD at the receivers. Note that reference fT4l where the authors considered Gaussian MIMO 
MACs with neither CSIT nor CDIT differs from our approach and that of HI, S, HOl, Ell, IM because in 
lfT4l the MIMO MAC is seen as a two-player zero-sum game where the first player is the group of transmitters 
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and the second player is the set of MIMO sub-channels. In the list of the aforementioned references, 191 seems 
to be the closest work to ours. However, our approach differs from |[9l on several technical key points. First of 
all, not only the BS but also the MSs can be equipped with multiple antennas. This is an important technical 
difference since the power control problem of jOl becomes a PA problem for which the precoding matrix of each 
user has to be determined. Also the issues regarding the existence and uniqueness of the network equilibrium 
are more complicated to be dealt with, as it will be seen. Specifically, random matrix theory will be exploited to 
determine the optimum eigenvalues of the precoding matrices. In Q, several assumptions made, especially the 
one involving the knowledge of all the instantaneous channels at each MS can be argued in some contexts. One 
of our objectives is to decrease the amount of signaling needed from the BS. This is why we assume that the 
BS can only send to the users sufficient training signals for them to know the statistics of the different channels 
and a simple and common coordination signal. The underlying coordination mechanism is simple because it 
consists in periodically sending the reahzation of a K!-state random signal, where K is the number of active 
users. As it will be seen in detail, such a mechanism is mandatory because, in contrast with ifTOll . we assume 
here successive interference cancellation (SIC) at the BS. Thus each user needs to know his decoding rank 
in order to adapt his PA policy to maximize the transmission rate. The coordination signal precisely indicates 
to all the users the decoding order employed by the receiver. Therefore the proposed formulation can be seen 
from two different standpoints. If the distribution of the coordination signal is fixed, then the addressed problem 
can be regarded as a non-cooperative game where the BS is imposed to follow the realizations of the random 
coordination signal. In this case the respective signal can be generated by any device (and not necessarily by 
the BS), in order to select the decoding order. On the other hand, if the distribution of the coordination signal 
can be optimized, the problem can be addressed as a Stackelberg game. Here the BS is the game leader and 
chooses his best mixed strategy (namely a distribution over the possible decoding orders) in order to maximize 
a certain utility, which will be chosen to be the network uplink sum-rate. 

In the described framework, one of our objectives is to know how well a non-cooperative but weakly 
coordinated system performs in terms of overall sum-rate w.r.t. its centralized counterpart (by "centralized" 
we mean that the users are imposed to follow the BS PA policies) when SIC is used at the BS. In this 
setting, several interesting questions arise. When the users' utility functions are chosen to be their individual 
transmission rates, is there a Nash equilibrium (NE) in the corresponding game and is it unique? What is the 
optimum way for a selfish user to allocate (spatially or temporally) his transmit power? How to choose the 
coordination signal that maximizes the network sum-rate? What is the performance loss of the decentralized 
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network w.r.t. the equivalent virtual MIMO network? 

This paper is structured as follows. After presenting the system model (Sec. Illl, we study in detail two PA 
games. In the first case (Sec. HID ), each MS is imposed to share his power uniformly between his transmit 
antennas but can freely allocate his power over time. In the second case (Sec. HVl ). we assume that the temporal 
PA is uniform and thus our objective is to derive the best spatial PA scheme. For each of these frameworks the 
existence, uniqueness, determination and sum-rate efficiency of the NE is investigated. Numerical results are 
provided in Sec. |V]to illustrate our theoretical analysis and in particular to better assess the sum-rate efficiency 
of the different games considered. We conclude the paper by several remarks and possible extensions of our 
work in Sec. |Vll 

II. System Model 

Throughout the paper v, M, (.)^ and (.)^ will stand for vector, matrix, transpose and transpose conjugate, 
respectively. For simplicity and without loss of generality, we will assume a MAC with K = 2 users. Note that 
the type of multiple access technique assumed corresponds to the one considered in the standard definition of 
the Gaussian MAC by |[T5l . |[T6l : all transmitters send at once and at different rates over the entire bandwidth. 
In this (information theoretic) context, very long codewords can be used and the receiver is not limited in 
terms of complexity. Thus the codewords of the different transmitters can be decoded jointly using a maximum 
likelihood decoding procedure (see fT6]| for more details). Interestingly, the transmission rates of the capacity 
region corresponding to the coding-decoding procedure just mentioned, can also be achieved, as discussed in 
lfT6l . by using perfect SIC at the receiver. In this paper we also adopt this decoding scheme, which means that 
not only the different channel matrices are perfectly known to the receiver but also that the codewords of all 
the users are decoded reliably. The case of imperfect CSIR and error propagation in the SIC procedure is thus 
seen as a useful extension of this paper. Since we assume SIC at the BS and that the users want to maximize 
their individual transmission rates, it is necessary for them to know the decoding order used by the BS. This 
is why we assume the existence of a source broadcasting a discrete coordination signal to all the terminals 
in presence. If this source is the BS itself, this induces a certain cost in terms of downlink signaling but the 
distribution of the coordination signal can then be optimized. On the other hand, if the coordination signal 
comes from an external source, e.g., an FM transmitter, the MSs can acquire their coordination signal for free 
in terms of downlink signaling. However this generally involves a certain sub-optimality in terms of uplink 
rate. Analyzing this kind of tradeoffs is precisely one of the goals of this paper. In both cases, the coordination 
signal will be represented by a Bemouilli random variable denoted with S ^ S. Since we study the 2— user 
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MAC, S = {1, 2} is a binary alphabet and S is distributed as Pr[5 = I] = p, Pr[5 = 2] = l-p = p. Without 
loss of generality we assume that when the realization of S is 1, user 1 is decoded in the second place and 
therefore sees no multiple access interference; in a real wireless system the frequency at which the realizations 
would be drawn is roughly proportional to the reciprocal of the channel coherence time (Tcoh)- Note that the 
proposed coordination mechanism is suboptimal in the sense that the coordination signal does not depend on 
the realizations of the channel matrices. We will see that the corresponding performance loss is in fact very 
small. 

We will further consider that each MS is equipped with rit antennas whereas the BS has antennas. In 
our analysis, the flat fading channel matrices of the different links vary from symbol vector to symbol vector. 
We assume that the receiver knows all the channel matrices whereas the transmitters have only access to the 
statistics of the different channels. At this point, the authors would like to re-emphasize their point of view: 

• On the one hand, we think that in some contexts our approach can be interesting in terms of signaling 
cost. We have seen that S lies in a element alphabet and the realizations are drawn approximatively 
at ^^[Hz], therefore the coordination mechanism requires at most bps from the BS and bps if 
it is built from an external source. Another source of signaling cost is the acquisition of the knowledge of 
the statistics of the uplink channels at the MSs. For example, in the context of coherent communications 
where the BS regularly sends some data to the MSs and channel reciprocity assumption is valid (e.g., in 
time division duplex systems) the corresponding cost can be reasonable. In general, this cost will have to 
be compared to the cost of the centralized system where the BS has to send accurate enough quantized 
versions of the (possibly large) precoding matrices at a certain frequency. 

• On the other hand, even if our approach is not interesting in terms of signaling, it can be very useful in 
contexts where terminals are autonomous and may have some selfish reasons to deviate from the centralized 
policies. In such scenarios, the concept of network equilibrium is of high importance. 

The equivalent baseband signal received by the BS can be written as: 



where xj^ (r) is the -dimensional column vector of symbols transmitted by user k at time r for the realization 
s € 5 of the coordination signal, Hfc(r) G c^^xnt channel matrix (stationary and ergodic process) of user 
k and z^*^(t) is an n^-dimensional complex white Gaussian noise distributed as M{0, (T^I„^); for sake of clarity 
we will omit the time index r from our notations. As ifTTl we assume that, for each s £ S, the data streams of 



K 




(1) 



k=l 
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user k are multiplexed in the eigen-directions of the matrix Q^,'*^ = E x)^'xX 

^ V^^^P^'^V^^^'^\ Finding 

(s) (s) 

the optimal eigen-values and coordinate systems V^, that maximize the transmission rate of user k is one 
of the main issues we will solve in the next two sections. In order to take into account the antenna correlation 
effects at the transmitters and receiver, we assume the different channel matrices to be structured according to 
the Kronecker propagation model ifTSl with common receive correlation |[T9l : 

yk G {1,...,K}, Hfe = Ri0feT| (2) 

where R is the receive antenna coiTclation matrix, Tj. is the transmit antenna correlation matrix for user k 
and @k is an x nt matrix whose entries are zero-mean independent and identically distributed complex 
Gaussian random variables with variance The motivation for assuming a channel model with common 
receive correlation is twofold. First, there exist some situations where this MIMO MAC model is realistic, the 
most simple situation being the case of no receive correlation i.e., R = I (see e.g., |[20l ). Although it is not 
explicitly stated in |[T9l the second feature of this model is that the overall channel matrix H = [Hi...Ha'] can 
also be factorized as a Kronecker model, which will allow us to re-exploit existing results from the random 
matrix theory literature. Therefore the case where the overall channel matrix is not separable can be seen as a 
possible extension of this paper that can be dealt with by using the results in 1(211 . 

In this paper we study in detail two special but useful cases of decentralized PA problems. In the first case 
(Game 1), we assume (for instance because of practical technical/complexity constraints) that each user is 
imposed to share his power uniformly between his transmit antennas but can freely allocate his power over 
time; this problem will be referred to as temporal PA game (Sec. Hill) . In the second case (Game 2), for every 
realization of the coordination signal, each user is assumed to transmit with the same total power (denoted by 
Pfc) but can freely share it between his antennas; this problem will be referred to as spatial PA game (Sec. ITVl ). 
For both games the strategy of user k € {1, 2} consists in choosing the distribution of for each s G 5 in 
order to maximize his utility function which is given by: 

s=l 

where 



EIog|I + r/HfeQ^^^Hf I if A: 



s 



(4) 



E log |I + r? ELi HfeQi'^Hf I - E log |I + r/H.^Q^tH^J if k ^ s 
with rj = and the usual notation for —k, which stands for the other user than k. Note that we implicitly 
assume Gaussian codebooks for the two users since this choice is optimum in terms of their individual Shannon 
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transmission rates (see e.g., 11221 ). This is why the strategy of a user boils down to choosing the best pair of 
covariance matrices (Q^.^"*, q[^^). The corresponding maximization is performed under the following transmit 
power constraint for each MS: Tr (^^^^^^ Pr[5' = s]Qi.^^) < UtPk- The main difference between Games 1 
and 2 relies precisely on how this general power constraint is specialized. In Game 1, the precoding matrices 
are imposed to have the following structure: V/c G {l,2},Vs G {1,2},q[,''^ = a^j^^Pklnt, which amounts to 
rewriting the total power constraint as follows 

2 

J]Pr[S = s]4'^ < 1. (5) 

s=l 

On the other hand, in Game 2, the power constraint expresses as 

VA: e {1, 2}, Vs e {1, 2}, Ti{q['^) < ntPk- (6) 

In both game frameworks, an important issue for a wireless network designer/owner is to know whether by 
leaving the users decide their PA by themselves, the network is going to operate at a given and predictable state. 
This precisely corresponds to the notion of a network equilibrium, a state from which no user has interest to 
deviate. The main issue is to know if there exists an equilibrium point, whether it is unique, how to determine 
the corresponding strategies and characterize the efficiency of this equilibrium in terms of network sum-rate. 

III. Temporal power allocation game 

As mentioned above, in the temporal power allocation (TPA) game, the strategy of user k G {1,2} merely 
consists in choosing the best pair (a[.^\a[^^). Since each transmission rate is a concave and non-decreasing 
function of the a[f'''s, each user will saturate the power constraint (5) i.e., X]^=i ^^I*^ ~ •^I'^i*^ ~ 1' which 
leads to optimizing a single parameter a^^^ or a^^\ From now on, for sake of clarity we will use the notations 
a^'^ = ai, = 02- Indeed, it is easy to verify that the power constraints are characterized completely, for 



the first user by = with ai G Aj^^ = 0, | , and for the second user by ag^^ = ^ (i^p)«2 ^^^^ 



(1) _ 



02 G Aj^^ = 0, . Thus the strategy of user k G {1,2} consists in choosing the best fraction ak from 
the action set -4^^^. Our main goal is to investigate if there exists an NE and determine the corresponding 
profile of strategies (a^^,a^^). It turns out that the issues of the existence and uniqueness of an NE can be 
properly dealt with by applying Theorems 1 and 2 of ||23l in our context. For making this paper sufficiently 
self-contained, we review here these two theorems (Theorem 2 is given for the 2— user case for simplicity and 
because it is sufficient under our assumptions). 

Theorem 1: lt23]| Let Q = {K,^{Ak}keic,{'Uk}keK) be a game where K, = {1^...,K} is the set of players, 
Ai, ■■■,Ak the corresponding sets of strategies and ui, Ufc the utilities of the different players. If the following 
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three conditions are satisfied: (i) each Uk is continuous in the vector of strategies (ai, ax) G ElfcLi -^fc/ (H) 
each Uk is concave in G .4^; (Hi) Ai, ■■■,Ak are compact and convex sets; then Q has at least one NE. 

Theorem 2: f23\ Consider the K -player concave game of Theorem\l\with K = 2. If the following (diagonally 
strict concavity) condition is met: for all (a'^, a") G A\ and {a^-, 0.2) € A^ such that (a'^, 7^ (o^i) (^2)^ Wi ~ 



> 0; then the uniqueness of the NE is 



m<^2) - £7K,«2)J +K-4) [ii^K,«'2) - i^K'«2) 

insured. 

At this point we can state the first main result of this paper, which is provided in the following theorem. For 
sake of clarity we will also use the notations: pk = p if k = 1 or pk = p if k = 2. 

Theorem 3 (Existence and uniqueness of an NE in Game 1): the temporal PA game described by: the set 
of players JC = {1,2}; the sets of actions A^^^^ = 0,^ and utilities u^i^ak, a_k) = pi?[,^^(afc, a_fc) + 
pR^^\ak,a-k)> where the rates R^^^ follow from Eq. (I?]) has a unique NE. 
Proof- 
Existence of an NE. It is guaranteed by the geometrical and topological properties of the utility functions 
and the strategy sets of the users (over which the maximization is performed). Indeed, we can apply |[23l 
in our matrix case. Without loss of generality, let us consider user 1. The utility of user 1 comprises two 
terms corresponding to the two coordination signal realizations: ui (01,02) = pR^i\-,-) +pRf'\-, ■). Using 
the fact that J^log|X + tYY^| = -Tr [Y" {X + tYY^y^YY^ {X + tYY^y^Y] it is easy to verify 
that ^^(01,02) = -ETr[BB^] < and ^^(01,02) = -ETr[CC^] < where B = piHf (I + 
piaiHiHf)-iHi, C = |piHf f I + pii^HiHf + /52a2H2H^) ~' Hi and pi = rjPi, p2 = VP2 corre- 



spond to the signal-to-noise ratios of the users. Thus for every user k, the utility is strictly concave w.r.t. to 
afc. Also it is continuous in (ai, over the convex and compact strategy sets Aj^^. Therefore the existence 
of at least one NE is guaranteed. Interestingly, we observe that for a fixed game rule, which is the value of 
the parameter p, there will always be an equilibrium. The users adapt their strategies to the rule of the game 
in order to optimize their individual transmission rates. 

Uniqueness of the NE. We always apply 1231 in our matrix case (see Appendix [All and prove that the 
diagonally strict concavity condition is actually met. The key of the proof is the following Lemma which is 
proven in Appendix |Bj 

Lemma 1: Let A', A", B' and B" be Hermitian and non-negative matrices such that either A' 7^ A" or 
B' 7^ B". Assume that the classical matrix order >z is total for each of the pairs of matrices (A', A") and 
(B', B") i.e., either A' t A" (resp. B' t B") or A" y A' (resp. B" y B'). Then we have Tr(M + N) > 
with M = (A" - A') [(I + A')~^ - (I + A")^^], N = (B" - B') [(I + B' + A')~^ - (I + B" + A")"^]- 
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It can be shown (see Appendix |A] for more details) that the diagonally strict concavity condition writes in our 
setup as pT(^) +pT(2) > where Vs G {1, 2}, T^*) is defined by T^**) = Tr(M(^) + N^*)) where the matrices 
M^''), N^'^) have exactly the same structure as M, N in the above Lemma. For example, if we consider two 
pairs of parameters {a'i,a'() £ [Aj^^)'^ and (02,02) G (^J^^^)^ such that either a'l ^ a" or a'g '^2 
in Theorem 121 T^^^ can be obtained by using the following matrices A' = pia'^HiH^^, A" = pia'/HiH{^, 



B' 



P2 



H2H|^, B" = p2 112^2 ■ The term T*^^^ has a similar form as T^^^ thus, applying Lemma 
[T] twice and considering the special structure of the four matrices (A', A", B', B"), one can prove that the 
term pT^^^ + pT^^'' is strictly positive. Therefore the unconditional uniqueness of the NE is guaranteed. ■ 
Determination of the NE. In order to determine the selfish PA of the users at the NE, we now exploit the 
large system approach derived in ll24l for single-user fading MIMO channels. This will lead us to simple 

approximations of the utility functions which are much easier to optimize. From now on, we assume the 

nt 

asymptotic regime in terms of the number of antennas: nt — > 00, — > 00, and lim — = c < 00. 

rif — >oo,nr— »oo n^. 

In this asymptotic regime, references |[24l . |[25l . |[26l provide an equivalent of the ergodic capacity of single-user 
MIMO channels, which corresponds exactly to the situation seen by user 1 (resp. 2) when 5 = 1 (resp. 5 = 2); 
this gives directly the approximation of the rates R^^^ and R^2^; see Eq. From Eq. dU we also see that 
the rates and R^2^ correspond to the difference between the sum-rate of the equivalent Knt x Ji^ virtual 



MIMO system and an nt x nr single-user MIMO system, therefore the results of 11241 . 11251 . 11261 can also be 
applied directly. The corresponding approximates can then be easily checked to be: 



R!{\ai,a2) 



-R2^^(ai, 02) 



1 + 7](iiPid\^>{i)^i I + ^log2 1 + Vd^^\j)5i 

1 - pa2 



i=l j=l 

nt _ nt 

J^log2 [l + 27?aiPidf)(i)72 

i=l 



i=l 



l + 2r/- 



P 



nj?77i(5i log2e 

(^)72 



(7) 



+ J^log2 [l + 2?yd(^)(j)52 



where V/c G {l,2},4^(i), i G {l,...,nt} are the eigenvalues of the transmit correlation matrices (see 
Eq. Q), d^^\j), j G {1, ...,71^} , are the eigenvalues of the receive correlation matrix R and the parameters 
7j, 6j are the unique solutions of the following systems of 2— degree equations: 



71 



61 



ntfr^^l + vd^){j)6i 
nt ^ 1 _L P. A^) I 



(8) 



^ l + r/aiPi<^(i)7i 
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( 1 " 

72 



1 ^ 



.2 ^ ^ 

2nt 



2nt ^ l + 27?d(R)(j),52 

1 + 27?aiPidp)(.)72 ^ 1 + 27?i^P24^^(i)72 



~ f 2^ ~ f 2) 

The approximate functions R\ (•,•) and i?2 (') can be obtained in a similar way and the approximated 
utility of user k G {1,2} follows: %(ai,a2) = pR^k\cii,Ci2) + 02) ■ Now, in order to solve the 

constrained optimization problem, we introduce the Lagrange multipliers (An, A12, A21, A22) G [0, +0x0)"^ and 
define for k G {1,2} the function ^^(ai, 02, A^i, Afc2) = -Uk{ai,a2) + A^i {au - ^) - Afc2afc- The Kuhn- 
Tucker optimality conditions follow. Therefore, the optimum selfish PAs, (a^^, a^^), can be obtained by using 
a fixed-point method and an iterative algorithm, following the same idea as in i fTOl for non-coordinated MIMO 
MACs with single-user decoding. At this point we have to make an important technical comment. Our proof 
for the existence and uniqueness of the NE holds for the exact game. For the approximated game, we need 
the approximated utilities to have the same properties as their exact counterparts. It turns out that the large 
system approximation of the ergodic mutual information can be shown to have the desired properties |[27l . In 
particular, the results of |[27l show that the approximated utilities are strictly concave and that if the iterative 
PA algorithm converges, it converges towards the global maximum. 

Sum-rate efficiency of the NE. Now, let us focus on the sum-rate of the decentralized network and compare it 
with the optimal sum-rate of its centralized counterpart. The centralized network sum-rate, denoted by Rsum, is 
by definition obtained by jointly maximizing the sum-rate over all the pairs of power fractions (ai, 02) € [0, 1]^: 
-^sum — loaax ui{ai,a2) +U2 (01,02). Knowing that log | • | is a concave function, one can easily verify 
that the maximum is obtained for (0^,02) = (1, 1) and that i?sum = Elog |l + piUiH^ + p2'H.2il2\- As 
the optimum precoding matrices are proportional to the identity matrix, it can be checked that the network 
sum-rate at the NE (denoted by i?^^) equal to the centralized network sum-rate for p = and p = 1: 
-^sum(O) = ^sJim(l) = -Rsum- Indeed, let us consider that p = 1. In this case, user 1 is always decoded in the 
second place (Pr[S' = 1] = 1). This means that there is no temporal power allocation game here and each user 
always allocates all of his available power for the case where 5 = 1: (a^^,a^^) = (1,0). Replacing in Eq. 
(4) the corresponding correlation matrices: q|^^^ = I^^, Q^^^ = I„j and Q^^-* = 0„j (the square zero matrix), 
= 0„, we obtain that R^^^{1) = R^L 

In the high SNR regime, where 77 — > cxo, we obtain from dUl,® that rjSi —>■ — , 77(^2 — > t— ^nd thus 71 and 



71' 272 

72 are the unique solutions of the following equations: ^ J2]Li ^.Td^^KJ) = 2^ Si=i -y2+d^%) = 
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approximated utilities become: 



rtt n,, 

B!i\ai,a2) = ^log2 1 + r/aiPi4'^^(i)7ij + ^log2 
i=i j=i 

R^^\ai,a2) = ^log2 [l + 2??aiPidp)(i)72] + ^log2 



1 + 



J J 



i=l 



1 + 27] 



71 
1 - pa2 



nt log2 e 



P 



1 + 



72 



-2nt log2e - i?^^^(ai,a2) 



(10) 



By setting the derivatives of •) w.r.t. ai and U2(-, •) w.r.t. 02 to zero, we obtain that, for each user, the 
PA at the NE is the uniform PA (aj^^, 02'^) = (1, 1)> regardless of the distribution of the coordination signal 
p G [0, 1]. Therefore, at the equilibrium, we have that 

-n-sumW — P^l l«l '0;2 j + P-K^ (a^ ,02 j+P-K2 ,02 )+PR2 i«l ' "2 J 

= pElog|I + piHiHf I +pElog|I + piHiHf + p2H2H^| -pElog|I + p2H2H|^| (11) 
+pElog |I + piHiHf + P2H2H|^| - pElog jl + piHiHf I + pElog |I + p2H2Hf | 



Knowing that the uniform spatial PA is optimal in the high SNR regime ifTTl . ifTOl . the centralized network 
sum-rate coincides with the sum-capacity of the centralized MAC channel, i?sum = C'sum- 

In the low SNR regime, where ^ 0, we obtain from dSjl, dill that 7]6i — > 0, 77152 —>■ and thus 71 = 
Y^]=i d^^Kj) ^ 72 = Y^]=i d^^\j) ■ Approximating ln(l + x) ^ x for x « 1, the achievable rates 
become: 

^ rir nt 



nt 



^2^''(ai, 02) 



1 1 
—Ti- 
nt p 



i=i 
pa2 



i=l 



nt 



(12) 



P2Y.S''Hj)Y.diii)^og2e 

j=i 1=1 

We see that the utilities Ukiai, 02) = ^vPk Sj=i d^ {i) log2 e converge and the network sum- 



nt"' 

rate at the NE coincides here again with the centralized network sum-rate: 



K 



(C) 



TTt d^^\j) [i]Pi Er=i 4 (0 + r]P2 Er=i 4 (^)j log2 e. In this case also, the price of anar- 

chy ll28l is minimal for any distribution of the coordination signal. 

To sum up we have seen that there is no loss of optimality in terms of sum-rate by decentralizing the PA 
procedure in at least four special cases: I) p = 0; 2) p = I; 3) when 7] co for any p G [0, 1]; 4) when r] 
for any p G [0, 1]. Additionally, in case 3), since there is no loss by imposing the spatially uniform PA lITTl . 
ifTOl . the centralized (and cooperative) MAC sum-capacity is achieved. If we further assume that there is no 
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correlation among the transmit antennas, = I, the uniform spatial PA is optimal lITTl for any rj. Thus, the 

(C) 

centralized sum-rate is always identical to the sum-capacity of the centralized MAC channel, i?sum = Csum- 
This means that if the BS chooses to use a completely unfair SIC-based decoding scheme, the selfish behavior 
of the users will always lead to achieving the centralized sum-capacity. This result is in agreement with O, 
where the authors have proposed a water-filling game for the fast fading SISO MAC (assuming perfect CSIT 
and CSIR) and shown that the equihbrium sum-rate is equal to the maximum sum-rate point of the capacity 
region. However, as opposed to the SISO MAC with the proposed coordination mechanism IH, the decentralized 
MIMO MAC with coordination does not achieve the sum-rate of the equivalent virtual MIMO network for any 
value of p and for an arbitrary noise level at the BS. In particular, the fair choice p = ^ is not optimal. We will 
quantify the corresponding performance gap through simulation results. Furthermore, in the low and high SNR 
regimes, the centralized sum-capacity is also achieved for any value of p. The consequence of these results is 
that any binary coordination signal can be used without loss of global optrmality. 



In this section, we assume that the users are free to share their transmit power between their antennas but 
for each realization of the coordination signal the transmit power is constrained by Eq. (6). In other words we 
assume that the users cannot distribute their power over time: they cannot decide the amount of power they 
dedicate to a given realization of the coordination signal. As a consequence of this power constraint (Eq. (6)), 
the two precoding matrices that each user needs to choose can be optimized independently and each of them 
does not depend on p. The strategy set of user k is in the spatial PA (SPA) game is: 



Theorem 4 (Existence and uniqueness of an NE in Game 2): The SPA game defined by the set of players 
K = {1,2}, the strategy sets A^^^^^ and utilities Uk{ak,a-k) given by Eq. (O, has a unique NE. 

Proof: The main feature of the game under the aforementioned power constraint is that there exists a unique 
NE in each sub-game defined by the realization of the coordination signal. The proof is much simpler than 
that of the time PA problem since the use of Rosen's Theorem |[23l is not required. Without loss of generality 
assume that 5 = 1. Whatever the strategy of user 2, user 1 sees no interference. Therefore he can choose 
q[^^ independently of user 2. Because r'^\q[^\ Q^^^ ) is a strictly concave function to be maximized over a 
convex set, there is a unique optimum strategy for user 1. As we assume a game with complete information 
and rational users, user 2 knows the utility of user 1 and thus the precoding matrix he will choose. The same 



IV. Spatial power allocation game 




(13) 
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concavity argument can be used for r'^\q^^\q'^^) and therefore guarantees that user 2 employs a unique 
precoding matrix. ■ 
Determination of the NE. In order to find the optimum covariance matrices, we proceed in the same way as 
described in ifTOl . First we focus on the optimum eigenvectors and then we determine the optimum eigenvalues 
by approximating the utility functions under the large system assumption. In order to determine the optimum 
eigenvectors, the proof in ll20l can be applied in our context to assert that there is no loss of optimality by 



is) (s) fj 

restricting the search for the optimum covariance matrix when imposing the structure = U^P^ XJj: , where 
Ufc is a unitary matrix coming from the spectral decomposition of transmit correlation matrix = Uf^DkUj^ 
defined in Eq. Q and the diagonal matrix P^j^^ = Dia.g{Plf\l), Plf\nt)) represents the powers user k 
allocates to the different eigenvectors. As a consequence, we can exploit once again the results of |[24l . ll25l . 
|[26l assuming the asymptotic regime in terms of the number of antennas. The new approximated rates are: 



i=l 
2 nt 



n,. 

J^l0g2 [l + 7?d(^)(i)5l 



J = l 



e=i i=i 



.(l)^.•^^(T), 



ntr/7i5i log2 e 

4nj??72(52 log2 e- 



(14) 

where VA; € {l,2},(i^ (z), i e {!,..., nj} are always the eigenvalues of the transmit correlation matrices T^, 
j € {l,...,nr} , are the eigenvalues of the receive correlation matrix R and the parameters ji, 5j are 
the unique solutions of the following systems of equations: 



71 
5i 



a)/,•^^(T)< 



(15) 



nt 



f,l + r,Pr{i)dY'{ihi 



72 



2nt ^ l + 27?(iW(j)(52 



2n 



1 ^ lit 

EE 



>(i) 



(16) 



^Sl + 27?p/^)«ci!^)«l72 



Then, optimizing the approximated rates R 



w.r.t. Pl^\i) leads to the following water-filling equations: 



VA; € {1,2}, 



{1),NE. 



ln2A 



(1) 



(17) 
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where a[.^^ >0, ke {1,2}, are the Lagrangian multipUers tuned in order to meet the power constraints given 

use the same iterative PA algorithm as the one described in ifTOl . Under 
the large systems assumption, in this game also, the approximated utilities have the same properties as the exact 
utilities. 

Sum-rate efficiency of the NE. Unlike the temporal PA game, we have not assumed a particular structure for 
the precoding matrices and thus the centralized solution coincides with the sum-capacity of the virtual MIMO 
network, iisum = C'sum- Another important point to notice here is that the equilibrium precoding matrices do 
not depend on p. This considerably simplifies the BS's choice for the sum-rate optimal value for p. Indeed, as 
we have already mentioned, the precoding matrices do no depend on p and therefore the sum-rate i?sum(p) is 
merely a linear function of p\ (p) = ap + b where 

a = Elog |I + r?HiQ«'^^Hf + r^H^Q^^^'^^H^I - Elog |I + r^HiQ? ^'^^Hf + r^H^Q^ '^^H^| 
b = Elog|I + r?HiQf)'^^Hf + r?H2Qf'^^H^|. 

(18) 

Depending on the sign of a, if the BS wants to maximize the sum-rate, it will choose either p = or p = I. 
If it wants a fair game it will choose p = ^ and accept a certain loss of global optimality. Note that even 
for p € {0, 1} the sum-capacity is not reached in general: this is because the matrix Q^^^'^^ (resp. Q2^^''^^) 
does not coincide with the first (resp. second) component of the pair of precoding matrices that maximizes the 
(strictly concave) network sum-rate. However, as we did for the temporal PA game, in the low and high SNR 
regimes one can show that the decentralized MIMO MAC has the same performance (w.r.t. the sum-rate) as 
its equivalent Knt x virtual MIMO network. 

V. SIMULATION EXAMPLES 

All the results will be provided by assuming the asymptotic regime in the numbers of antennas. We know, 
from many contributions (see e.g., |[TOl . ||27l . |[29l . ||30l ) that large-system approximates of ergodic rates are 
accurate even for relatively small systems. We also assume that R = I. 

For the TPA problem, we look at the case where there is no transmit coiTclation, T^. = I. We have seen that 
the performance of decentralized MAC depends on the rule of the game i.e., the value of p. This is exactly 
what Fig. [T] depicts for the following scenario: Pi = 1, P2 = 10, r? = 5 dB, = = 4. First, we see that 
the MAC sum-rate is a convex function of p and the maximum of -R^m(p) is reached for p G {0, 1}. In these 
points, which correspond to the most unfair decoding schemes (either user 1 or 2 is always decoded first) the 
centraUzed sum-capacity of the MAC is achieved. One important observation to be made is that the minimum 
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and maximum only differ by about 1%. Many other simulations have confirmed this observation. This shows 
that whatever the value of p, the gap between the sum-rate of a decentralized MIMO MAC with selfish users 
and the sum-capacity of the equivalent cooperative MAC (virtual MIMO network) is in fact very small. Now, 
we want to evaluate the benefits brought by using a SIC instead of single-user decoding lITOl . For the scenario 
where Pi = P, P2 = lOP with P G [0, 20], = nj = 4 and r/ = 5 dB, Fig. 2 shows the achievable network 
sum-rate at the NE versus the available power at the first transmitter P. For the SUD scheme, the users are 
decoded simultaneously at the receiver. In this case both users see all the interference coming from the others. 
We see that the SIC scheme performs much better than the proposed SUD scheme, regardless of the distribution 
of the coordination signal: this comparison makes sense especially for the point p = ^ since both decoding 
schemes are fair. 

From now on, we consider the SPA problem. In this case we assume an exponential correlation profile for 
Tfc such that Tk{i,j) = t^^ (note that Tr(Tfc) = nt), where < tk < I is the corresponding correlation 
coefficient |[3ll . Il32l . We already know that the sum-rate is a linear function of p and therefore is maximized 
when either p = or p = 1. It turns out that this slope has a small value. Furthermore, it has been observed 
to be even for a symmetric MAC, i.e.. Pi = P2 and ti = t2. These observations have been confirmed by 
many simulations. In Fig. 3 we have plotted the sum-rate achieved by varying p for the scenario: Pi = 5, 
P2 = 50, r/ = 3 dB, rit = Ur = 4, ti = 0.4, t2 = 0.3. Even in this scenario, which was thought to be a 
bad case in terms of sub-optimality, the sum-rate is not far from the sum-capacity of the centralized MAC. 
For the same scenario, we have plotted in Fig. 4 the achievable rate region and compared it to that obtained 
with SUD. We observe that in large MIMO MAC channels, the capacity region comprises a full cooperation 
segment (approximately) just like SISO MAC channels. The coordination signal allows one to move along an 
almost straight line, corresponding to a relatively large range of rates. 

VI. CONCLUSION 

We have provided complete proofs for the existence and uniqueness of an NE in fast fading MIMO MACs 
with CSIR and CDIT where the transmission rate is chosen as user utility. By exploiting random matrix theory, 
we have also provided the corresponding optimum selfish PA policies. We have seen that the BS can, through a 
single parameter (i.e., p G [0, 1], which represents the distribution of the coordination signal), force the system 
to operate at many different points that correspond to a relatively large range of achievable transmission rate 
pairs. We know, from [H, 191 that for Gaussian MACs with single antenna terminals, this set of rate pairs 
corresponds to the full cooperation segment of the centralized MAC. Said otherwise a decentralized Gaussian 
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SISO MAC with coordination achieves the same rate pairs as a MAC with full cooperation or virtual MIMO 
system. The goal here was to know to what extent this key result is valid for fading MAC with multi-antenna 
terminals. It turns out this is almost true in the MIMO setting. In the cases of interest, where the power is 
optimally allocated either over space or time, the performance gap is relatively small even though the proposed 
coordination mechanism was a priori sub-optimal since it does take into account the channel realizations (known 
to the receiver). Interestingly in large MIMO MACs, the capacity region comprises a full cooperation segment 
just like SISO MACs. The coordination signal precisely allows one to move along the corresponding (almost) 
straight line. This shows the relevance of large systems in decentralized networks since they allow to determine 
the capacity region of certain systems whereas it is unknown in the finite setting. Furthermore, they induce an 
averaging effect, which makes the users' behavior predictable. Indeed, in large MIMO MACs the knowledge of 
the CSIT does not improve the performance w.r.t. the case with CDIT. To conclude we review some extensions 
of this work which we have suggested throughout it. It would be interesting to study the case of the decentralized 
space-time PA, which, in particular, would require the generalization of Lemma [T] to arbitrary positive matrices 
and exploitation of some results in |[2T1- A second useful extension would be to evaluate the impact of a non- 
perfect SIC on the PA problem. At last, we will mention that it would be useful to evaluate analytically or 
bounding the price of anarchy of the NE, which would require to find a bounding technique different from that 
used for non-atomic games |[34l . |[35l . |[36l . 



Appendix A 

We want to prove that the diagonally strict concavity condition is met for the time PA problem i.e., for all 

{a[,ai) £ [Aj^-^y and (02,02) G [Aj^'^)'^ such that either o[ 7^ a" or 02 7^ 02 we want to prove that: 



C = {o'l-o[) 



dRi , , , , dRi . „ „^ 

("1,02) - TT— («li"2. 



doi ' doi 
We can write C = pT^^^ + pT^"^^ where for all s G {1, 2}: 



Or 



dR2 , , ,s 



> 0. 



T(^) = (a'/ -a'l) 



do\ 

By expanding T*^^) we have 



{0^,02) - -^^[0^,02) 



+ («2 - «2) 



'-'^2 I II II \ 



(19) 



(20) 
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TW = K-a;)ETr{[(I + pia;HiHf)-i-(I + pia;HiHf)->iHiHf} 

+{a'i - a'2)ETV {(I + pialHiHf + l^p^Yl^il^)-^ p^^n^Hf^ 

- (I + pi«'/HiHf + lzMp2H2H^)"V2f H2H|^} 

= ETr {(A" - A')[(I + A')"^ - (I + A")"^] + (B" - B')[(I + B' + A')"^ - (I + B" + A")"^]} , 

(21) 

where A' = pia'iHiHf , A" = pia'/HiHf , B' = pz^^f^HaH^, B" = pa^^f^HaH^. We observe that 
the matrices A', A", B' and B" verify the assumptions of Lemma[T] First, they are Hermitian and non-negative. 
Second, as they write as A' = a'HiHf , A" = a"HiHf , B' = ft'HaH^ and B" = 6"H2H^, we also see 
that the matrix order ^ is total for each of the pairs of matrices (A', A") and (B',B"). This directly follows 
from the fact that the scalar order > is total, which implies that either a" > a' or a" < a' and either b" > b' or 
b" < b'. By considering the particular structure of the four matrices and applying Lemma [T] it is straightforward 
to see that the term T^^^ is strictly positive, T^^) > 0. In a similar way we can prove that T^^) > and thus 
the diagonally strict concavity condition is met: C > 0. 

Appendix B 

Proving Lemma [T] amounts to showing that 

T = Tr {(A - B)(B-i - A^^) + (C - D)[(B + D)~i - (A + C)-^]} > (22) 

where the matrices A = I + A", B = I + A', C = B" and D = B' have been introduced for more clarity. 
Since the matrix order >z is total for A and B, and C and D it suffices to prove that T > for the four 
following cases: (1) A ^ B and C ^ D; (2) A ^ B and C ^ D; (3) A ^ B and C ^ D; (4) A ^ B and 
C ^ D. 

Case (1): A ^ B and C ^ D. To prove the desired result in this case we use the following lemma. 

Lemma 2: If ^ is a Hermitian and non-negative f M = ^ j and N is non-negative (N ^ j but not 
necessarily Hermitian, then Tr(MN) > 0. 

Proof: We write Tr(MN) = Tr(M^/2NM^/2) > 0. We have used the fact that M is a Hermitian non- 
negative matrix to write M = M^/^M^/^. Knowing that N is a non-negative matrix one can easily check 
that M^/^NM^/^ is also a non-negative matrix and thus the trace (sum of the non-negative eigenvalues) is 
non-negative. ■ 
The quantity T writes as T = Tr(MiNi) + Tr(M2N2) where Mi = A-B, Ni = B'^ - A-\ M2 = C-D 
and N2 = (B + D)^^ — (A + C)~^. Clearly these four matrices are Hermitian. Since by assumption Mi >z 
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and M2 ^ we only need to verify that Ni >z and N2 ^ to be able to apply Lemma |2] to T. The matrix 
Ni is non-negative because for any pair of invertible matrices (X, Y): X ^ Y -<4> Y^^ >z X^^ (see e.g., |[33il ). 
The same result applies to N2 since by assumption A + C ^ B + D. Using lemma |2] concludes the proof. 

Case (3): A ^ B and C ^ D. To treat this case we first prove the following auxiliary Lemma. 

Lemma 3: Let X and Y be two distinct, Hermitian and positive matrices of size n: X = X^ >- 0, Y = 
yO andX^ Y. Then Tr[(X - Y)(Y-i - X'^)] > 0. 

Proof: It is easy to see that Tr[(X - Y)(Y"i - X^^)] = Tr[Z + Z^^ - 21], with the Hermitian and 



n 



positive matrix Z = Xi Y^^Xi and thus we further have Tr[(X - Y)(Y-i - X^^)] = V ^^^^^\ ^' > 

where the matrix Az = Diag(Az(l), Az(n)) corresponds to the spectral decomposition of . ■ 
By applying this lemma to T we have that: 

T = Tr{(A-B)(B-i-A-i) + [(C + A)-(B + D)][(B + D)-i-(A + C)-i] 

-(A-B)[(B + D)-i-(A + C)-i]} (23) 
> TV {(A - B)(B-i - A"i) - (A - B)[(B + - (A + C)-^]} . 

We know that C ^ D then C + A ^ D + A and also that (C + A)~^ ^ (D + A)~^. Using the fact that 
A ^ B and also Lemma |2] we have that Tr[(A - B)(C + A)~^] > Tr[(A - B)(D + A)~i] and the trace 
becomes lower bounded as T > Tr {(A - B)(B~i - A~^) - (A - B)[(B + D)"^ - (A + D)"^]}. Now, we 
are going to prove that this lower bound, say 7/,^, is positive: 

Tlb = Tr{(A-B)(B-i-A-i)-[(A + D)-(B + D)][(B + D)-i-(A + D)-i]} ^^^^ 
= TV |(A - B)(B-i - A"i) - [(A + I) - (B + I)] [(B + I)-^ - (A + } 

where we have made the following change of variables: A = D^/^AD^/^, B = D^/^BD^/^ such that A = 
D-1/2AD-1/2 = A^ ^ and B = D-i/^BD-^s = ^ q. By applying the Woodbury formula (A + 
= A~i - A-^{A + 1)^1 and (B + 1)"^ = B^^ - B^{B + the lower bound T^^^^ rewrites as: 



\2 



Tlb = TV|(A-B) 



B-i - A-i - B-i + B-ifB + 1)^1 + A^i - A^ifA + I) 



TV |aB-1(B + + BA-i(A + - (A + - (B + . 



(25) 



Let us denote the ordered eigenvalues of the two matrices A and B as A^(l) < A^(2) < . . . < A^(n) and 
A^(l) < A^(2) < . . . < A^(n). From ||37l we know that for two matrices X and Y of size n, Tr(XY) > 
Er=i Ax(«)Ay(ra - i + 1), which implies directly that Tr(XY'i) > Y17=i where Ax(«) and Ay(f) are 

the ordered eigenvalues (in the previously specified order) of the corresponding matrices. Applying this result 
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we find that 

Tr [aB"1(I + B)-i 

and finally obtain that: 



> 



-A^(i)(l + A^(i))^ 



TV 



BA-i(I + A)-i 



> 



E 



> ^ [A^W-A^(i)]^[l + A^(i) + A^(i)] 



(26) 



(27) 



A^(i)Aj^(i)[l + A^(i)][l + A5(i)] 
To conclude the global proof one can easily check that Case (2) (resp. Case (4)) can be readily proved from 
the proof of Case (1) (resp. Case (3)) by interchanging the role of A and B and C and D. 
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Fig. 1. Temporal PA game. Achievable network sum-rate versus p for Pi = 1, P2 ~ 10, n,- — rit = 4, i] — 5 dB. The sum-capacity 
of fading MIMO MACs is reached for both unfair SIC decoding schemes (pi = and P2 = 1) and is very close to this upper bound 
for any distribution of the coordination signal, Vp G (0, 1). 
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Fig. 2. Temporal PA game. MAC sum-rate versus the transmit power P for Pi — P, P2 = lOP, = nt = 4, r; = 5 dB. Comparison 
between the fair SIC decoding scheme (p — i), the unfair SIC scheme (p = 0), and SUD decoding scheme. 
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Fig. 3. Spatial PA game. MAC sum-rate versus p for Pi = 5, P2 = 50, rir — nt = i, rj = 3 dB, ti — 0.4, t2 = 0.3. The achievable 
network sum-rate of fading MIMO MACs is linear w.r.t. p G [0, 1] and is very close to the centralized upper bound. The optimal 
distribution obtained with the Stackelberg game is p* = 0. 
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Fig. 4. Spatial PA game. Achievable rate region for Pi = 5, P2 = 50, Ur = nt = 4, rj — 3 dB, ti — 0.4, t2 — 0.3. By varying p 
allows to move along a segment close to the centralized sum-capacity, similar to the SISO MAC channels. 
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